NASA/TM— 2000-210032 



Sensitivity of a Wave Structure 
to Initial Conditions 

Walter M.B. Duval 

Glenn Research Center, Cleveland, Ohio 


May 2000 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STT Program Office provides access to the 
NASA STT Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STT Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 

has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
data bases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STT 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http:! fwww.stiMasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at (301) 621-0134 

• Telephone the NASA Access Help Desk at 
(301) 621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 


NASA/TM— 2000-210032 



Sensitivity of a Wave Structure 
to Initial Conditions 

Walter M.B. Duval 

Glenn Research Center, Cleveland, Ohio 


Prepared for the 

38th Aerospace Sciences Meeting and Exhibit 

sponsored by the American Institute of Aeronautics and Astronautics 
Reno, Nevada, January 11-14, 2000 


National Aeronautics and 
Space Administration 


Glenn Research Center 


May 2000 



Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 
Price Code: A03 


National Technical Information Sendee 
5285 Port Royal Road 
Springfield, VA 22100 
Price Code: A03 


SENSITIVITY OF A WAVE STRUCTURE TO INITIAL CONDITIONS 


Walter M.B. Duval 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

Microgravity experiments aimed at quantifying 
effects of g-jitter via controlled sinusoidal forcing 
transmitted on the interface between two miscible 
liquids have shown the evolution of a quasi -stationary 
four-mode wave structure oriented vertically. The sensi- 
tivity of the wave structure to phase angle variation is 
investigated computationally. We show that a slight 
variation of the phase angle is sufficient to cause a 
bifurcation to a two-mode structure. The dependence of 
phase angle on wave structure is attributed to sensitivity 
on initial conditions due to the strong nonlinearity of the 
coupled field equations for the parametric space of 
interest. 

Introduction 

Transport processes which are generated by buoy- 
ancy induced flows due to nonconservative body forces 
occur through an interface. In many practical problems 
such as crystal growth of binary 7 alloys via directional 
solidification 1 or solution crystal growth 2 the distribu- 
tion of the concentration field near the interface deter- 
mines the quality of the resulting crystal. This class of 
problem involves the coupling of the density field, 
through the body force, to the concentration and or tem- 
perature field. The availability of a low gravity environ- 
ment, such as a Space Shuttle or the International Space 
Station, offers a laboratory 7 in w r hich the body force is 
reduced to such an extent that mass transport occurring 
at an interface can potentially occur uniformly. How- 
ever, random vibration or g-jitter in that environment 
can affect mass transport so as to negate the advantage 
of the low r body force. We study a model problem 
which considers the effect of controlled body forces to 
lead insight in understanding effects of g-jitter. 

We consider the problem of two miscible liquids 
juxtaposed inside a cavity, with a sharp density interface 
initially, such that the body force is parallel to the inter- 
face. Experimental study to understand the effect of g- 
jitter on such a system, using the Microgravity Vibration 
Isolation Mount (MIM) to input controlled vibration 
forces, 3 shows that the interface can become unstable 
against Kelvin-Helmholtz instability and generate quasi- 


stationary waves with various modes. The effect of con- 
trolled vibration input to the system is to produce an 
oscillating shear flow at the interface. This problem has 
similarity to the classical unbounded parallel shear flow 
4.5.6.7.8.9.10.11 j n s tably stratified fluids which generate 
Kelvin-Helmholtz instability. We investigate computa- 
tionally the sensitivity of quasi-stationary (q-s) waves on 
orientation of the oscillatory component of the body 
force through phase angle variation for the parametric 
space in which these structures were observed. We 
show that the observed two-mode q-s waves are 
independent of phase angle for low Stokes-Reynolds 
( Re s) number. However, when the effect of nonlinearity 
becomes important for the four-mode q-s waves, there 
exists narrow regions of phase angle for which the 
four-mode is stable. 

In the following w ; e investigate computationally 
the stability of two and four mode q-s waves when the 
initial condition on the phase angle is varied. We show 
the dynamics of the flow r field which produces the q-s 
waves as well as the magnitude of the velocity field as a 
function of phase angle. For the parametric region in 
which nonlinearity becomes important we deduce the 
range of phase angle for which the four-mode q-s w^ave 
is stable. 

Model Problem 

Consider two miscible liquids introduced side by 
side inside a cavity as shown in Figure la. Such 
configuration is possible under a low r gravity environ- 
ment since the hydrostatic pressure jump can be made 
negligibly small due to the low steady component of the 
body force n g a . For certain parametric regions the 
interface can be stable and remain stationary in the vert- 
ical orientation. The effect of g-jitter may be studied 
systematically by introducing controlled vibration inputs 
rng o cos((O/+0) parallel to the interface. This can permit 
the study of the effect of amplitude (m,g 0 ), frequency 
(Co) and phase angle (0) on the system. For the 
parametric range (ng 0 ) for which the interface is stable, 
one may investigate vibration inputs that cause 
instability. 
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Given a stable interface oriented vertically inside a 
cavity separating two miscible liquids, the system will 
deviate from thermodynamic equilibrium; the interface 
will spread over time such that the system becomes 
homogeneous or uniformly mixed over a long time scale 
on the order of (L 2 /D AB) . However depending on the 
body force field g(t), for short time scales, the interface 
can remain stable or form various patterns according to 
the species continuity equation 

nr 

~ = D AB v-c. (1) 

The initial value of the two liquids may be prescribed as 
follows, 


must be satisfied. The overbar denotes average values. 
The body force considered consists of a steady and 
oscillator} 7 component of the form 


f(t) = n g a 4- mg a cos(cor + 9) 


( 10 ) 


t - 0 


C(.v,y,0) = 


1 0 < x < L/l 

.5 x = L/l (2) 

0 L/l <x < L 


This class of problem represents an initial value prob- 
lem, whose solution may depend sensitively on the ini- 
tial conditions set by the prescribed density of the fluids 
or the nonconservative body force field prescribed by 
n g 0 , mg a , and 0; 0 may also be varied independently 
while keeping the other variables fixed to investigate the 
sensitivity to initial conditions. 

Dimensionless parameters of the problem may be 
obtained by selecting the characteristic length, time, 
velocity, as 


- T = 

‘ V CO ’ 


I . „ Ap mgoLc 2 

— and U c ~ — - 

CO p v 


The fluid on the left and right hand side has respectively 
the value of 1 and 0, and the interface has the value of 
0.5, The boundaries have a normal (/?) and the condition 
of impermeability is satisfied, 


VC rt =0 on T. 


( 3 ) 


For parametric space in which the nonlinear convective 
terms in equation (1) are negligible, if the cavity domain 
is approximated to be infinite, a closed form solution 
shows 




(4) 


The Stokes length scale was selected since the growth of 
waves is generated by the short wavelength instabilities, 
12 and it also leads to a minimum set of dimensionless 
numbers to characterize the problem. When the steady 
component of the body force dominates, the geometric 
length scale such as the height (H) and length (L) 
becomes more important. 

For computation the above set of equations are 
transformed to vorticity-stream function formulation by 
taking the curl of the momentum equation, let 

dV -d*¥ ^ = dv du 


u — 


v = 


f AB l 


that the concentration field diffuses uniformly as a func- 
tion of time. Whereas for parametric regions in which 
the interface becomes unstable, the nonlinear convective 
terms become dominant. 1 
teristic of the flow field 

\? =?(x,yj- A) (5) 

dictates the pattern of the interface for a vector space of 
parameters A. Since the density field, which is a linear 
function of the concentration, field 

p = p(l + (3 AC) (6) 

is coupled to the body force, the momentum equation 

= -Vp + pV 2 \? + psV) (7) 


dy ' dx 

The resulting set of equations become 

d 2 ^ a 2 *? 

a^ + 


dx dy 


(II) 


means that the charac- 

4- Re, 

a; ac 


d% 

dt 

dx dy 


. dx 2 dy 2 . 


( AR + cos( t + 0 ))-^ 

ox 


ac 

3/ 


+ Re v 


dC dC 


i 

5c 


d 2 C d 2 C 
+ 


a.v 2 dv 2 


( 12 ) 
■ G3) 


Since we are considering a cavity with aspect ratio 1 
(H=L), the parametric space can be succintly represented 
as 


A = A(Re s , AR , Sc , 0) 


and the equation of incompressibility 

vV = o 

with no slip at the boundaries (V) 

i? = o on r 


in which, 


( 8 ) 


( 9 ) 


Re v= 


Ap m So 

s=— v> — 1/2 

p C0 v “V 


II c V 

, AR = — , Sc =— — 

III D 
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Res denotes a Stokes-Reynolds number, amplitude ratio 
AR, and Schmidt number Sc. 

Equations 11-13 were discretized in space using 
finite difference techniques. The sharp discontinuity of 
the interface was resolved using the Flux Corrected 
Transport algorithm for computation. We used third 
order Adams-Bashforth for time discretization. All com- 
putations were carried out for a 90x90 grid size. We 
are interested in investigating the sensitivity of phase 
angle variation on the solution of the above equations 
for a parametric range of interest. 

Discussion and Computation Results 

Parametric Range 

We consider the stability of two and four-mode 
quasi -stationary waves as a function of phase angle for 
two cases. The two-mode case corresponds to viscous 
fluid pairs consisting of 22% deuterium oxide (D?(9) 
with de-ionized water (HiO) on the left hand side, and 
100% de-ionized water on the right side. The kinematic 
viscosities are equal, v=0.01 cm 2 /sec, with a binary dif- 
fusion coefficient D AB =2xlO~ 5 cm 2 /sec, and density ratio 
A(yjT= 0.022. The body force has a steady component 
n ^=lxl0" 6 (go , Earth’s gravitational acceleration), and 
oscillatory component m= 1.72x10^ with a frequency 
f=0.1 Hz. The values of the dimensionless parameters 
are Re s =2.4xl0'\ AR= 5.8xl0 _ \ Sc=500 (l/Sc=2xlCT'). 
The four-mode case consists of a dilute mixture of 
ethylene glycol (20%) and water on the left side and 
100% de-ionized water on the right, with a mean viscos- 
ity v=0.0108cm 2 /s, D AB =lx\0^an 2 /s , and A(yjT=0.022. 
The steady component of the body force is the same as 
above, however, m=2xl0 -2 or 20m g and f=l Hz. The 
dimensionless parameters are respectively, Re 5 =0.3375, 
AR=5x 10'\ Sc=1079 (1/Sc=9.3xl0~ 4 ). The cavity size 
is of height and width H=L=5 cm, hence aspect ratio of 
1. The phase angle was varied from 0 to 27t for each 
case. 

Global Bifurcation of the Interface 

We first discuss the base state case for which the 
phase angle is zero. Under ideal microgravity condi- 
tions in which all disturbances are isolated, it is possible 
to obtain an initial condition as shown in Figure la. 
The interface band is represented by values of .1 to .9 
so that the effect of mass diffusion can be shown as it 
evolves. Since the two fluids are miscible, the interface 
would spread over time while remaining stationary in 
the vertical orientation to satisfy the species continuity 
equation (1) for Re*=0 over a time scale T~L~/D AB of 
order 0(5x1 0 4 sec) for a 1 cm diffusion width. How- 
ever, when energy is injected into the system through 
some controlled means, such as a MIM, the characteris- 


tic time for instability to grow is on the order of 100 
sec, depending on the parametric input. For excitation 
with a frequency of 0.1 Hz, the interface becomes 
unstable against Kelvin- Helmholtz instability due to an 
oscillating parallel shear flow which results in a two- 
mode quasi-stationary wave, Figure lb. The oscillating 
shear flow is generated by the oscillatory component of 
the body force. The instability is largely due to the 
creation of vorticity at the interface as shown in the last 
term of equation (12). Note that since the amplitude 
ratio is small AR=0.0058, the mean field is negligible as 
a result the interface remains stationary, otherwise it 
would tilt. The spreading or width W of the isoconcen- 
trates which indicates mass diffusion equal to W=0.55 
cm is due vibration driven convection. For comparison, 
the width predicted for mass diffusion when convection 
is negligible from equation (4) 

5 = 2 Vd as / (14) 

is on the order of 0.28 cm which is approximately 0.5 
times less than the convective case. The nonlinear con- 
vective terms in the concentration field equation (13) 
nearly balance the diffusion terms 
(Re s =0.0024= VSc =0.002), thus one would expect mass 
diffusion to occur. 

Increasing the effect of nonlinearity (ftc 5 =0.3375) 
as seen from the vorticity equation (12), increases the 
number of modes to four, Figure lc. The intensity of 
the flow field drives the instability. In this case the 
interface remains relatively sharp in comparison to the 
two-mode case. The concentration field equation (13) 
shows that the nonlinear convective terms are dominant 
over the mass diffusion term 

(7?c s =O.3375:»l/Sc=9.3xl0~ 4 ). This means that the inter- 
face can be deformed by the flow field for the time scale 
considered with minimum mass diffusion. The interface 
width W=0.22 cm is 1.58 times greater than the diffu- 
sion distance 5=0.126 cm. 

Flow Field Dynamics as a function of Phase Angle 

Two -mode Case 

The flow field dynamics, which yields the two- 
mode structure, is shown in Figure 2. Various phase 
angles are shown which typifie trends for variation in 
the range of 0° to 360". As the phase angle is varied 
the wavelength 

X = X(Re Sf AR , 5c, 0) 

of the q-s wave remain constant, independent of phase 
angle. The jump discontinuity of the density at the 
interface initially causes the interface to evolve like a 
vortex sheet. The tangential discontinuity of the velo- 
city field causes the oblate shape of the streamlines. 
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Diffusion of mass at the interface smears the discon- 
tinuity as the flow evolves. The two-mode structure is 
produced by a single vortex. When the phase angle 
changes to 90° the effect of the steady component of the 
body force dominates which effectively changes the flow 
field to a rotational flow. Figure 2b. This can be seen 
from the vorticity equation, note that at 90° the oscilla- 
tory component of the body force is zero initially, thus 
the steady component dominates. This example under- 
scores the sensitivity to initial condition of the flow 
field. On the other hand, for a phase angle of 180° 
which corresponds to the condition when the oscillatory 
component of the body force goes through a relative 
minimum, Figure 2c, the single cell structure bifurcates 
to multi-cells. The flow field is weak, as indicated by 
the magnitude of the velocity field, thus there is no 
effect on the deformation of the interface. Even though 
the flow field is unstable with respect to phase angle 
variation, the interface deformation remains stable and 
independent of phase angle. 

The magnitude of the velocity field as a function 
of phase angle is shown in Figure 3. The behavior of 
the flow field is not symmetric with respect to 180, 
there is a decrease in the magnitude of the flow field for 
the phase angle range of 180" -360". The maximum 
magnitude of the velocity occurs at 90° and 270" w'hile 
the minimum occurs at 0" , 180", and 360°. This exam- 
ple illustrates that for low Re y (0.0024), even though 
there is a 200% variation in the magnitude of the velo- 
city field, the wavelength of the interface remains con- 
stant independent of the phase angle. 

Dynamics of Flow Field for four-mode Structure 

The dynamics of the flow field that generates a 
four-mode quasi -stationary wave is shown in Figure 4 
for a phase angle of zero degrees. Initially there is an 
oscillating parallel shear flow at the interface, Figure 4a, 
the tangential discontinuity of the velocity field is shown 
in the vector field plot. Unlike the unbounded parallel 
shear flows 5 which generates travelling Kelvin- 
Helmholtz instability waves, the boundary in our case 
cancels out the wave motion and produces quasi - 
stationary 7 wave patterns. The recirculation of the flow r 
field due to impingement and reflection at the top and 
bottom boundaries create a region of concentrated vorti- 
city, thus secondary cells are formed. The interface 
evolves like a vortex sheet due to discontinuity of the 
tangential velocity at the interface region, Figure 4b. A 
line vortex is created due to Kelvin-Helmholtz instability 
that results in undulation of the interface. The vortices 
grow as time increase, Figure 4c, two row's are formed 
similar to unbounded shear flows as shown by the vector 
field. The crests and troughs of the interface is due to 
co-rotating vortex pairs, as showm by the vector field, 


which forms homoclinic orbits. The singularity point 
between the vortex pairs serve as natural centers which 
stretch and rotate neighborhood portions of the interface 
into its undulatory shape. 

Stability of the four-mode Structure as a function of 
Phase Angle 

The effect of phase angle on the wavelength of 
the interface as well as the dynamics of the flow is 
showm in Figure 5. The basic four-mode structure as 
observed experimentally'* occurs for a phase angle of 0" , 
175°, and 360". Typical bifurcation of the four-mode 
structure is shown for 9 = 4°, 90", and 270". Note 
that a slight change of phase angle of 4" causes a transi- 
tion to a two- mode structure. For phase angle in the 
range of 5"-172" the initial characteristic of the flow 
field changes from an oscillating parallel shear to a rota- 
tional flow with a single cell which has a stable elliptic 
point. However, the oscillator}’ component of the body 
force destabilizes the elliptic point to a homoclinic 
point. Thus, two co-rotating vortices are formed as 
shown by the streamlines and vector field in Figure 5c 
for 90°, The signature of the oscillatory parallel shear 
flow' is preserved. The co-rotating vortices effectively 
stretch and deform the interface into its sinusoidal pat- 
tern. A clockwise rotation of the flow field results for 
phase angle in the range of 5° -172", Whereas, for 
phase angle in the range of 180" -355" a counter- 
clockwise rotation of the flow occurs, Figure 5d. 

The magnitude of the velocity field as a function 
of phase angle is showm in Figure 6. Unlike the two 
mode case in Figure 3, there is symmetry with respect 
to 180". The maximum velocity occurs at 90" and 270", 
whereas the minimum occurs in the neighborhood of 
173" and 353° respectively. The basic four-mode struc- 
ture occurs within 3" bands, for phase angle of 0"-3", 
175" -178", and 357" -360". The regions of stability of 
the four-mode structure is represented by the average 
dimensionless wavelength X*=(A/7/) as a function of 
phase angle in Figure 7. The narrow' dip of the curve 
where the relative minimum occurs represents regions of 
decreasing wavelength or higher mode numbers. 
Overall the four-mode structure is sensitive to phase 
angle, it is stable over a narrow’ neighborhood region, 
near 0", 175", and 360". Within the bands of 0"-5", 
175" -180", 355" -360", shorter wavelengths occur. The 
trends show that for increasing Stokes-Reynolds number 
Res, the wavelength of the interface depends on phase 
angle. The effect of nonlinearity is to shift the 
wavelength of the interface to shorter and shorter 
wavelengths. The results show the importance of initial 
conditions on wave structure and the importance of 
aligning an experiment correctly with respect to the 
body force in a microgravity environment. 


NASA/TM — 2000-210032 


4 



Conclusions 

The stability of two and four-mode quasi- 
stationary waves is investigated as a function of initial 
conditions for phase angle in the range of 0° to 360°. 
Two-mode q-s waves which are formed at low Stokes- 
Reynolds numbers (Re^ =0.0024) are stable with respect 
to phase angle variation even though the flow field is 
unstable; the wavelength is independent of phase angle. 
However, when nonlinearity becomes important for 
higher Re 5 (0.3375), the four-mode wave bifurcates to a 
two-mode over most of the phase angle range. Only a 
narrow region exists for which the four-mode is stable. 
The wavelength is dependent on phase angle. This 
finding underscores the sensitivity to initial conditions of 
experiments performed for parametric spaces where non- 
linearity becomes important. 
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